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(54) Method and apparatus for reducing three-dimensional shape data 



(57) A method and apparatus for reducing three-di- 
mensional shape data, includes calculating estimation 
values for portions to be deformed by converging two or 
more vertices (Vj) of a polygon model (M,) f and reducing 
the number of data for the polygon model (Mj) by con- 
verging two or more vertices (Vj) of the polygon model 
(Mj) when the estimation values are equal to or below 
the predetermined permissible value. Estimation values 



for surfaces to be deformed by shrinking edges or sur- 
faces are calculated based on distances between the 
respective surfaces and all the original vertices (v,) in- 
volved in the surface deformation. Before each data re- 
duction, the portion having been involved in all the pre- 
vious data reductions is defined as a reduction prohibi- 
tion area, and a succeeding data reduction is applied to 
a portion other than the reduction prohibition area. 
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Description 

[0 0011 Thisinventionrelatestoam^ 
SS A huge number of measurement data are obtained 

'apparatus having quite a high space resolution if an .magequ *. h ^S^Lng«io«^|Koo«^ 
quality than necessary, itdisadvantageousiy results,^ ^£^£22 this proble m as much as possible by 
and a considerable reduction in process.ng E^JTTJiTSr?. drffe rence between a polygon model (here- 
reducing the number of measurement data. Howev e rj Inn su* a data after the data reduction and a 
inafter, referred to as approximated polygon mod °? (hereinafter, referred to as original 
polygon mode, formed by the ™ asured J^ cannot be effective* utilized. 

xr^n^ 

SoTpaSaV 
STC^b^ 

using the actually measured three-dimensiona,^ "^^ESuo. the number of three-dimensional data 
[00041 On the other hand, vanous methods have ^ ^ IdTwS a desired permissible error range, 
forming a triangular polygon model by approxtmat ^^^S^Smi Polygon Models" (by Horikawa and Tot- 
[0005] For example, ™r^£$^^ ^ disc.os^a method for reducing the number 

are short. . innathan Amitabh Varshney, Dinesh Manocha, Greg Turk, 

[00061 Further, in "Simplification Envelopes" (by Cohen Jonathan ™" aDn Proce edings of SIGGRAPH 96 in 
Hans Weber. Pankai Agarwal. Frederick Brooks ^^^^^S^d^.tn^terg^ 
Computer Graphics Proceedings. Annual Conferen^ Sen^ 

Carnegie Mellon Unh/ersrty, Proceedings of S '^^ H J;'_^ by converging two or more vertices (three-d.men- 
35 three-dimensional shape data fbrnmng a x " an ^°*J°™ one vertex. In this method, a distance 

sional shape data) forming an edge or surface of ^ 

between theconverged vertex due to edge sh^^^^ as an error ln the data 

.o g^ 

Kot'h^^ 

45 dimensional shape data within a P r f at *^^ of Three-Dimensional Polygon Mod- 

[0009] Sincethedatareducingmethoddsclosedin H ™'^3 the ^^ation of a model by the edge 

els" oniy deletes edges based on evaluation 
deletioncannotbe estimated although thethree-dimens^ 

Thus, this method is effective in the case of ^ u ™.^ this method cannot be 

error. However, it is necessary to generate polygon B^ 1 ^ 9 number of data relating to 

STSSir. ft* mefto* dMd In -Su«* «*■— Us,„ 9 Qu.** a. Metrics?" conft* ft. <M-~ 
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from the vertices of the polygon model after the approximation from the surfaces forming the original polygon model 
as errors. Thus, this method is not necessarily a suitable method in the case that data reduction is desired to be 
performed by controlling errors of the polygon models after the approximation from the vertices forming the original 
polygon model like the measurement data of the three-dimensional shape measuring apparatus. 
5 [0012] Further, in the case of measurement data obtained by measuring rough surfaces or cut surfaces of a meas- 
urement object by means of a three-dimensional shape measuring apparatus, e.g., in the case that the orientation of 
surfaces locally suddenly change in a zigzag manner in a polygon model, errors for these zigzag surfaces are estimated. 
Therefore, errors are estimated as large values as the data reduction progresses, thereby presenting the problem of 
insufficient data reduction. 

10 [0013] Furthermore, each time the edge or surface shrinking is performed, it is necessary to estimate errors with 

respect to the edges or surfaces of the shrunk polygon model that are to be influenced by the edge or surface shrinking. 

This will increase the number of calculations, consequently causing the problem of hindering the data reduction. 

[0014] It is an object of the present invention to provide a method and apparatus for reducing three-dimensional 

shape data which are free from the problems residing in the prior art. 
15 [0015] It is another object of the present invention to provide a three-dimensional shape data reducing method and 

apparatus which enable higher speed processing. 

[0016] It is another object of the present invention is to provide a three-dimensional shape data reducing method 
and apparatus which enable proper data reduction based on precise error estimation for measurement data obtained 
by a three-dimensional shape measuring apparatus. 
20 [0017] According io an aspect of the invention, a method for reducing three-dimensional shape data, comprises 
calculating estimation values for portions to be deformed by converging two or more vertices of a polygon model, and 
reducing the number of data for the polygon model by converging two or more vertices of the polygon model when the 
estimation values are equal to or below the predetermined permissible value. 

[0018] The portions to be deformed may be surfaces of the polygon model. Estimation values for surfaces to be 
25 deformed by shrinking edges or surfaces may be calculated based on distances between the respective surfaces and 
all the original vertices involved in the surface deformation. 

[0019] Alternatively, the number of data for the polygon model may be reduced by converging two or more vertices 
of one portion of the polygon model after another portion, repeatedly. Before each data reduction, the portion having 
been involved in all the previous data reductions is defined as a reduction prohibition area, and a succeeding data 

30 reduction is applied to a portion other than the reduction prohibition area. 

[0020] According to another aspect of the invention, an apparatus for reducing three-dimensional shape data, com- 
prises an estimation value calculator for calculating estimation values for portions to be deformed by converging two 
or more vertices of a polygon model, and a data reducing device for reducing the number of data for the polygon model 
by converging two or more vertices of the polygon model when the estimation values are equal to or below the prede- 

35 termined permissible value. 

[0021] The estimation value calculator may calculate estimation values for surfaces to be deformed by shrinking 
edges or surfaces based on distances between the respective surfaces and all the original vertices involved in the 
surface deformation. 

[0022] Also, the data reducing device may reduce the number of data for the polygon model by converging two or 
40 more vertices of one portion of the polygon model after another portion, repeatedly. Further, there may be provided a 
judging device forjudging based on a calculated estimated value whether a portion of the polygon model is permissible 
for the deformation, and a prohibition area defining device for defining a portion having been involved in all the previous 
data reductions executed by the data reducing device as a prohibition area to keep a succeeding data reduction from 
being applied. 

45 [0023] These and other objects, features and advantages of the present invention will become more apparent upon 
a reading of the following detailed description and accompanying drawings, in which: 

FIG. 1 is a block diagram showing a construction of a three-dimensional shape measurement data processing 
system according to an embodiment of the invention; 
so FIG. 2 is a block diagram showing a construction of a shape data processor; 

FIG. 3 is a diagram showing an example of approximation of a polygon model by data reduction; 
FIG. 4 is a concept diagram showing a data reducing method by edge shrinking; 
FIG. 5 is a concept diagram showing a data reducing method by surface shrinking; 

FIG. 6 is a diagram showing an estimation value used in a data reducing method according to the embodiment; 
55 FIG. 7 is a flowchart showing a procedure of a three-dimensional shape data reduction processing; 

FIG. 8 is a diagram showing edge shrinking in a portion where a discontinuous surface is formed on a polygon 
model by a random noise data; 

FIG. 9 is a diagram showing a difference between an error estimation value according to the embodiment and a 
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HGMS is a diagram showing a modification of the flowchart shown in FIG. 13. 

menslonal shape data of a measurement object G to be measured ^ ^erna. g proC essor 3 

dimensional shape data obtained by the thjee-dimensional shape MJP^ Mishap i ^ 

for reducing the three-dimensional shape data obtained by the 4 for storing the three-di- 

data reducing method to perform a specified data Processing an ^^ ^^^^Z.^^ 
mensionalshape data processed by shape data processor 3 ^ after data 

image of the measurement object G such as a po^gon moden u»ng ^^SSSi and computer aided 
reduction, and an external processing dev.ce 6 such as a ^^~ m P^ 8 ^ data p r0 cessor3 areutilized. 

=eT b ;rr u rn^ 

G on a measurement table 1 06 provided in front of the ^"^^ the measurement 

the surface shape of the measurement object G at eac^ 

table 106 at speeded intervals along height directs. ^S^SS!m^L^^ a direction normal 
height direction of the measurement object G, the opfcal S^asSnt object G is measured by 

to £ direction and Z-direction, respectK,e.y, the ^^^%^^^oZ^Q tne Y . C oordinate of the 
measuring Z-coordinate of the measuremenl ^bjectG ^^^^ ^Zwessu^t^G 
measurement object G at specified .ntervals. The three- dimens £ * ' 3 
measured by the three-dimensional shape ?J£3S5L* £5K*«* 'converting device 

[0027] The three-dimensional shape measuring apparatus 1 .s P^iaeow"" p > rf jtion 

03a arranged in parallel to X-axis and formed of a line sensor and the ^JSSwS^lSS. corresponds 
of the measurement object G is measured using a l^ 

of the slit extending in X-direction is emitted from the '»™' n ^^ 0 f the optica, system 

of the optical system 102. This illuminat.on ght s , ga *ered at an objec ^L^dy^ plM-!0 oiectrtc€my coiivertlnQ 
102 to formaline beam extending in X-direct.cn at this ^8P^»J^^'™P , ne ^ )rttte 
device 103a is repeated* exposed in a sp eclf .c cycle wh, j h e 22J*2J3^^^ 0(lh . nM „ 

of the measurement object G. . „ Q ~w*i a data reducina device 32, a permissible 

[0029] Theshapedataprocessor3isprovidedw.thapo^ 
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three-dimensional shape measuring apparatus 1 . The data reducing device 32 reduces the three-dimensional shape 
data by edge shrinking or surface shrinking of the polygon model. In other words, as shown in FIG. 3, a precise polygon 
model A based on the measurement data is approximated to a polygon model B within a specified permissible error 
range by suitably reducing the vertices (three-dimensional shape data) of the polygon. 

5 [0030] The data reducing device 32 includes an estimation value calculator 321 , a estimation value sorter 322 and 
an edge/surface shrinking section 323 to carry out a data reduction processing to be described later. 
[0031] The estimation value calculator 321 calculates estimation values used to judge whether the respective edges 
or surfaces can be shrunk before shrinking edges or surfaces forming the polygon model by converging two or more 
vertices to one vertex. This estimation value is described later. 

10 [0032] The estimation value sorter 322 includes a calculator for calculating a minimum estimation value by rearrang- 
ing a plurality of estimation values calculated for the respective edges or surfaces by the estimation value calculator 
321 in an increasing or decreasing order. The edge/surface shrinking section 323 compares the minimum estimation 
value calculated by the calculator of the estimation value sorter 322 with a permissible value inputted via the permissible 
value input device 33 and performs an edge or surface shrinking corresponding to the minimum estimation value when 

15 the minimum estimation value is equal to or below the permissible value. The estimation values and an edge shrinking 
method are described later. 

[0033] The permissible value input device 33 is operated by an operator to input the permissible value used to judge 
whether the edges or surfaces can be shrunk as described above. This permissible value corresponds to a permissible 
range within which the edge or surface shrinking is assumed to have a little influence on characterizing portion(s) of 
20 the polygon model in approximating the polygon model, and is suitably inputted by the operator according to the shape 
of the measurement object G. 

[0034] The data output device 34 is an interface for outputting the reduced three-dimensional shape data (three- 
dimensional shape data forming the approximated polygon data) to the external storage device 4, the display device 
5 and the external processing device 6. 

25 [0035] At feast part of devices of the shape data processor 3 may be provided in another independent processing> 
unit, and their functions may be realized by cooperation of the processing units and the shape data processor 3 
Alternatively, these devices may be processing units which are separately provided in the shape data processor 3 and 
cooperate with each other. Conversely, at least two functions of these devices may be realized by a processing program 
^of one processor, e.g., the shape data processor 3 or a processor provided therein. 

30 [0036] The external storage device 4 is adapted to store the data in an external storage medium such as a magnetic 
disk, an optical disk or a magnetooptical disk. Upon input of the three-dimensional shape data from the shape data 
processor 3, the external storage device 4 stores the data in an external storage medium in a specified file format. 
■[0037] The display device 5 is an electronic display device such as a cathode-ray tube or a liquid crystal display 
device. Upon input of the three-dimensional shape data from the shape data processor 3, an image of the approximated 

35 polygon model of the measurement object G is displayed on the display device 5 based on the inputted data (see the 
polygon model B of FIG. 3). 

[0038] The external processing device 6 is a device for fabricating, for example, a mold or a model having substantially 
the same configuration as the surface configuration of the measurement object G. Upon input of the three-dimensional 
shape data from the shape data processor 3, a mold or the like is automatically fabricated by cutting a work based on 
40 the inputted three-dimensional shape data. 

[0039] Referring back to FIG. 1 , at least part of the measuring apparatus 1 and the devices 2, 4, 5, 6 may be incor- 
porated in the shape data processor 3. 

[0040] Referring to FIGS. 4 and 5, next, the data reducing method will be described. The data reducing method by 
edge shrinking is such a method for reducing the number of data that, in the case that a polygon model is defined by 
45 vertices v 1 to v 10 as shown in FIG. 4, two vertices v n , v 2 are, for example, so converged to a vertex v 1 I as to shrink an 
edge Ed 12 connecting the vertices v 1t v 2 . 

[0041 ] I n the data reducing method by edge shrinking (hereinafter, referred to as "edge shrinking method"), estimation 
values for estimating the influence of the edge shrinking on a change of the configuration of the polygon model for all 
edges Ed g (i, j = 1, 2, ...10, combinations of vertices adjacent to each other when h*j) are first calculated, and the 

50 vertices v,, Vj at the opposite ends of the edge Edy corresponding to a minimum estimation value are converged to a 
specified vertex v{ for data reduction if the minimum estimation value is equal to or below a permissible value. 
[0042] Subsequently, the edge shrinking is performed in a similar manner for the polygon model (approximated 
polygon model) in which the vertices v jf Vj are converted to the vertex vf, thereby further reducing the data. Thereafter, 
the edge shrinking is performed edge by edge in a similar method until no more edge shrinking can be performed for 

55 the polygon model after edge shrinking (i.e., the approximation of the polygon model is repeated until the minimum 
estimation value becomes larger than the permissible value) to reduce the data. 

[0043] It should be noted that the vertex v,' is an approximated point set in a specified position between the vertices 
v,, Vj such as a middle point of the edge Edy as disclosed in the aforementioned paper "Hierarchical Approximation of 
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Three-Dimensional Polygon Models". 

[0044] The data reducing method by surface shrinking (hereinafter, referred to as merely "surface shrinking method ) 
is such a method for reducing the number of data that, in the case that a polygon model is defined by vertices v 1 to 
v„ as shown in FIG. 5, vertices v 1 , v 2 , v 3 are so converged to a vertex v,' as to delete a surface P 123 enclosed by the 
5 vertices v 1 , v 2 , v 3 . 

[0045] Basically, the data reducing method by surface shrinking is carried out in a procedure similar to that of the 
data reducing method by edge shrinking. Specifically, estimation values for estimating the influence of the surface 
shrinking on a change of the configuration of the polygon model for all surfaces P, jk (i, j, k = 1 , 2, ...10, combinations 
of vertices adjacent to each other when i*j*k) are first calculated, and the vertices v,, v j( v k defining the surface P ijk 
10 corresponding to a minimum estimation value are converged to a specified vertex v,' for data reduction if the minimum 
estimation value is equal to or below a permissible value. Thereafter, the surface shrinking is performed surface by 
surface in a similar method until no more surface shrinking can be performed for the polygon model after surface 
shrinking to reduce the data. 

[0046] In the case that the three-dimensional shape data are successively reduced by the edge reducing method or 

is surface reducing method, the three-dimensional shape daia forming the measurement object G after the data reduction 
processing are different from those of the measurement object G actually measured by the three-dimensional shape 
measuring apparatus 1 and serve as approximated three-dimensional shape data. Thus, it is necessary to perform the 
data reduction processing while controlling an error of the polygon model after the data reduction (hereinafter, referred 
to as "approximated polygon model") from the polygon model formed by the actually measured three-dimensional 

20 shape data (hereinafter, referred to as "original polygon model"). 

[0047] In the data reducing method according to this embodiment, if it is assumed that a tnangular polygon P^ is 
approximated to a triangular polygon P^' by, for example, approximating vertices a, b, c (actually measured three- 
dimensional shape data) to vertices a', b\ C as shown in FIG. 6 through, for example, repetition of the edge shnnking, 
distances of a surface including the triangular polygon P abc ' after the approximation from the original vertices a, b, c 

25 involved in the approximation of this surface are controlled as errors. 

[0048] In the example of FIG. 6, the vertices involved in the approximation are assumed to be the vertices a, b, c in 
order to simplify the description. However, in an actual processing, vertices other than the vertex a are involved in the 
approximation, for example, while the vertex a is converged to the vertex a'. Therefore, the above errors include dis- 
tances of these vertices from the triangular polygon P^'. 

30 [0049] Next, the estimation value of the error used in the data reducing method according to this embodiment is 
described. Generally, a surface p in a three-dimensional space is expressed by Equation (1): 

A-x+By+C-z+D = 0 0) 

where 



35 



40 



45 



50 



55 



[0050] On the other hand, if (v x , vy, v z ) denotes the coordinates of point v in the three-dimensional space, a distance 
d between the point v and a surface p is calculated by Equation (2): 

d=A - v x + B- v y + C v z + D (2) 
and can also be expressed by matrices as in Equation (3): 

d=P T V (3) 

where 
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a 
b 
c 
d 



10 pT is a transposed matrix of P. 

[0051] Accordingly, the square of the distance d is expressed by Equation (4): 

<? = (P T . V) 2 = (P T - V)iV T - P) = P T (V- V T )' P = P T • Kv • P 



15 



(4) 



where lA is a transposed matrix of V 
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Kv = V*V T 
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Y x Y y 
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v x v y 
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Y y K x 



V V 
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[0052] Equation (4) expresses the square of the distance d between one plane p and one point v not included in the 
plane p. In the case that a plurality of such points V 1 , v 2 , ...v m exist for the plane p, a sum A(p) of the squares of these 
points V| (i = 1 , 2, ...m) from the plane p is expressed by Equation (5): 



35 



•Kv r P) = P 7 -gjSVP 



...(5) 



40 



[0053] In the data reducing method according to this embodiment, if there are n planes of a polygon model which 
are influenced by the edge shrinking, for example, when data are reduced by edge shrinking, distances of a plurality 
of original vertices (original three-dimensional shape data) involved in the data reduction processing for the approxi- 
mation from the planes are estimated for all such planes. Accordingly, the sum A(p) of the squares of the distances is 
calculated by Equation (5) for each of such planes, and e which is obtained by dividing a sum of the sums A (p,) by a 
total number M of all vertices involved in the data reduction processing as seen in Equation (6) is assumed as an error 
estimation value: 



45 



50 



M 



...(6) 



55 



[0054] Further, Equation (7) is obtained by substituting 



7 



EP 1 113 400 A2 



of Equation (5) for A(P,) in Equation (6): 




It should be noted that in Equation (7) , 

P, denotes a column vector for an equation of a plane p, subject to the deformation; 

SCn^eT™ 

data of an original vertex v involved in the deformation of the plane Pl is defined by Kv = V ■ V( Vdenotes a column 

^denltefa'stm of error matrix Kv y for a vertex Vj involved in only the deformation of the i-th plane p,; and 
n denotes a sum of planes p, to be deformed. 

20 [0055] In Equation (6), A(p,) is a sum of the squares of all original vertices involved in the approximation of the i-th 
plane p, from the plane p„ and m, is the number of terms in the sum of the squares of the dist « n ^ s _ 
r00561 In the calculation by Equation (6), since a matrix Kv, is a symmetric matrix as shown m Equation (4) only 

gl elements and elemeVabove them may be stored in the memory for 
vertex v, This enables a memory capacity necessary for the calculation to be reduced. ^S^SI^^ 

25 m , in the denominator of Equation (6) are the number of terms of the sum «?i?ZZZ^VZ^2i£i 
approximated plane p.. Since the element at the right bottom comer of the matnx Kv is constently 1 and the element 
auhe right bottom comer of a matrix XKv in the ca.cu.ation of the sum of the squares of the 

by Equation (6) may be conducted using the value of this element. Thus, the calculation is easier because rt is not 

necessary to separately calculate the number of terms m,. _ ^.u, th .„ 

30 [00571 Furthe ^ it may be appreciated to use weighted error matrix for the calculation of the ermre.Specrf^ 

are provided a weight Ltor w, for each vertex V, in consideration of the influence of the vertex V, to an >PP»*n£d 
polygon model or other causes. The weighted error matrix is obtained by multiplying the error matrix Kv, for the vertex 
V, by a predetermined weight factor w, for the vertex V, as w, ■ Kv,. „ imntKloll 

[00581 lnthecasethatEquation(5)isappliedine^^ 

3S n the processing of edge shrinking or surface shrinking for the plane p are specified, a sum ZKv, of the .matrices Kv, 
obtained by the three-dimensional shape data of these vertices v, serves as a parameter mdicating a degree of the 
error in the edge shrinking or surface shrinking. In the description below. IKv, is rtf erred to , - , n 
r00591 Next a procedure of the three-dimensional shape data reduction processing is descnbed with reference to 
FIG 7 The tnree-Smensional shape data of the measurement object G obtained by the three-dimensional shape 

40 inputted^ the shape data processor 3, in which a triangular polygon -ode. ^ene^ted 

l^tS pZon generator 31 (Step #1). Further, an operator inputs a permissible value used to judge whether the edge 
shrinking processing can be performed via the permissible value input device 33 (Step #3) 
[0060] Subsequently, estimation values ej (j = 1 , 2, ...N, N is a total number of edges) ,n the case th at a. I edges of 
the orUal polygon model are shrunk are calculated by Equation (6) in the estimation value ca culator 321 1 (Step #5) 

45 For instance XL example of FIG. 4. estimation values ^ = 1,2, 19) in the case o shnnking *eedg« i Ed^r- 
, tn c o in! Ed is - 5 6 9) and Ed-.* (t = 3, 4, ...10, where t+1 = 11 is assumed to be t+1 =3) are calcuiatea 
rooeii At his « the edge Ed 12 is^mnk by. for example, converging the vertices v„ v 2 to the vertex v V . the 
fSSSL P^a'e ai. dSormS andlight pianes p, p 2 , ... Pl0 (hereinafter, p.anes after the ^^J^S 
to as approximated planes) are formed by the vertex v,' and the vertices v 3 to v 10 not involved in the edge sMnkmg 

so As3mesquaresofthedistancesA(p0(k=1.2....8)are 

Tota number M of the elements included in all the sums A (Pk ) Is calculated, ^ Stan (? 

shrinking of the edge Ed 12 fe calculated by Equation (6). The evaluat.on value e(Ed 12 ) is expressed by Equat.on (8). 



55 



{a( P ,)+A(j> 2 ) + " + A(p,)} 
v nJ (m, + m 2 +-+/n 8 ) 
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[0062] At this time, the plane p 1 is defined by the vertices v,', v 3 , v 4 . Since the vertices involved in the edge shrinking 
are only v 1f v 2 if the vertices v 1t v 2> v 3 , v 4 are actually measured data, the sum A^) of the squares of the distances 
is calculated by Equation (5) using an error matrix (Kv^Kvg) relating to these vertices v-, , v 2 . If, in this state, the vertex 
v,' and another vertex v m (unillustrated actually measured data) are converged to a vertex v^ (not shown) to approx- 

5 imate the approximated plane p t further to an approximated plane p t ' (not shown), the vertices involved in the edge 
shrinking leading to the approximated plane p t " are v 1f v 2 , v m . Thus, a sum A(p v ) of the squares of the distances is 
calculated by Equation (5) using an error matrix (Kv-j+Kv^KvJ relating to these vertices v 1l v 2 , v m . 
[0063] In the second edge shrinking, it is assumed that an error matrix Kv n ' for the vertex v^ is (Kvj+Kv^, and the 
sum A(p n ') of the squares of the distances is calculated by Equation (5) using the error matrix Kv^ for the vertex v.,' 

10 and an error matrix Kv m for the vertex v m . 

[0064] Accordingly, a general theory of the edge shrinking processing is as follows. If the vertices v 1} v 2 , v 3 , v 4 are 
data generated by the edge shrinking instead of being actually measured data in FIG. 4, the error matrix Kv t ' for the 
vertex v^ defining the approximated plane p t is Kv t '= Kv^+Kvj, when the error matrices for the respective vertices v 1t 
v 2 , v 3 , v 4 are assumed to be Kv, , Kv 2 , Kv 3 , Kv 4 . The sum A(p0 of the squares of the distances Is calculated by Equation 

is (5) using an error matrix Kp.,(= Kv 1 '+Kv 3 +Kv 4 = Kv^Kvg+KVg+Kv^ obtained by adding the error matrices Kv^, Kv 3 , 
Kv 4 for the three vertex v^, v 3> v 4 defining the approximated plane p v 
[0065] An equation used to calculate the sum A(p.,) of the squares of the distances is: 

2 o A(p 1 )=P r r .K p1 .P 1 (9) 

Here, P t is a column vector of the approximated plane p 1 and m 1 is a total number of the error matrices Kv included 
in the error matrix Kp^ 

[0066] For the approximated planes p 2 , p 3 , ... p 8 , sums A(p 2 ), Afp^, ... A(p 8 ) of the squares of the distances are 
25 calculated by the same method as for the aforementioned approximated plane p 1f and are substituted into Equation 
(8) to calculate the error estimation value e(Ed 12 ) in the shrinking of the edge Ed 12 . 

[0067] The error estimation values E(Ed 1r ), e(Ed 2s ), e(Ed tt+1 ) in the case that the other edges Ed 1r , Ed 2s , Ed tt+1 are 
shrunk by the same method are calculated. 

[0068] Subsequently, the error estimation values ej are sorted in the estimation value sorter 322 (Step #7) and it is 
30 discriminated whether the minimum estimation value Min{ej} is equal to or below the permissible value inputted in Step 
#3 (Step #9). If the minimum estimation value Min{ej} exceeds the permissible value (NO in Step #9), the data reduction 
processing by the edge shrinking is completed since a change of the shape resulting from application of the edge 
shrinking to all the edges is impermissibly large. On the other hand, if the minimum estimation value Minfej} is equal 
to or below the permissible value (YES in Step #9), the edge corresponding to the minimum estimation value Minfcj} 
35 is shrunk in the edge/surface shrinking section 323 (Step #11). For example, if the minimum estimation value Min{e]} 
is e(Ed 12 ) in FIG. 4, the edge Ed 12 is shrunk as shown in FIG. 4. 

[0069] Subsequently, this routine returns to Step #7 after error estimation values e q (q = 1 , 2, .. .) are calculated for 
the edges of the polygon model after the edge shrinking that have been influenced by the edge shrinking (Step #1 3). 
A minimum estimation value M\r\{e^ of these estimation values e q is calculated and the edge corresponding thereto is 
40 shrunk (Steps #7 to #11). Thereafter, the edge shrinking is repeated edge by edge (loop of Steps #7 to #13), and the 
data reduction processing is completed when the minimum estimation value Minfe^ exceeds the permissible value 
(NO in Step #9). 

[0070] By performing the data reduction processing as above, the polygon model after the data reduction is approx- 
imated from the original polygon model such that the error estimation values from the original measurement data are 

45 equal to or below the permissible value. 

[0071] As described above, in the data reduction method according to this embodiment, the distances of the planes 
of the polygon model approximated by the edge shrinking or surface shrinking from the original vertices (measured 
three-dimensional shape data) are treated as the error estimation values. Thus, the data reduction processing can be 
performed by the estimation standard intuitively easily understandable. 

so [0072] Further, since the data reduction is performed such that the evaluation values lie within the specified permis- 
sible error range, it can be performed at a high speed while relatively precisely controlling the errors of the three- 
dimensional shape data forming the polygon model after the approximation. 

[0073] Particularly, reduction in the reliability of the error evaluation values can be suppressed even in the case that 
random noise data are mixed into the measured three-dimensional shape data or data reduction is performed by 
55 converging a multitude of vertices to one point. 

[0074] FIGS. 8 and 9 show a comparison of the error estimating method according to this embodiment and the 
conventional error estimating method in the case that random noise data are mixed into the measured three-dimen- 
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JETS .^hows the edge shrinking of a portion (shown i, , se^on, where 

on a polygon model by random noise data. If vertxes v, v, of the ^porf on where large ^ ^ 

the outer surface o, the pofrgon ^™£Z^*?£^ ^ nlding vertices v,. v 4 and a plane p 25 
based on, for example, distances d, , , d 2 verte J * £ ^ ucing P m 4 e thoddisclosedintheabovepaper"Surface 

includin g verticesv 2 ,v 5 asshown.nFIG.9accort.ng to»e^n«uang estimation values are smaller than sensory 

latter. . ^ „„. oet imatina method according to this embodiment and 

[0076] FIGS. 10, 11 and 12 show a comparison of ™™?»^J^ m conve ^ed to one point, 
the conventional error estimating method in the v " v * " v « of 3 

rOOTT] FIG. 10 shows the shrinking of a z.gzag surface (sh own J^ nllw in the middle bet ween 
polygon model. If the vertices v,. v 2 . v, n J ^,^^15^^ -<V of the vertex v V 

the vertices v 7 and v 8 as shown in FIG. 10, error estimation , mode| M shown in FIG. 12 

to the respective zigzag planes p 71 , p 12 . Pa* P*s< ^ . Surface Simplification Using Quadric Error Met- 
according to the data reducing method disclosed m ^^^^^^ errors. Contrary to this, ac- 

the estimation values e for the edges influenced by the edge shrinking cha 
sM n«n g ™o=L,ng « b„ pert.™.,. „a the pen^ «. £pp, «£ £ W « . ^ - ^ 

rirp= 9 9 »r:t^ 

is prohibited (Step #1 3a) and then Step #7 follows shrinki na processing is performed while excluding 

performed (loop of Steps #7, #9, #11, #13a). 
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[0087] On the other hand, if the minimum estimation value Min{8q} becomes larger than the permissible value (NO 
in Step #9) while the edge shrinking processing is repeated, it is discriminated whether the edge shrinking has been 
prohibited for these edges (Step #1 5). If the discrimination result is affirmative (YES in Step #1 5) , Step #7 follows after 
error estimation values e, (r = 1 , 2, ...) are recalculated for these edges (Step #1 7). A minimum estimation value Min 
5 {e,} of the error estimation values are calculated, and the edge corresponding to the minimum estimation value Min 
{e r } is shrunk (Steps #7, #9, #11). 

[0088] If the minimum estimation value Minfeq} becomes larger than the permissible value and there is no more edge 
for which the edge shrinking is prohibited (NO in Steps #9, #15), the data reduction processing by the edge shrinking 
is completed. 

10 [0089] FIG. 14 shows an exemplary data reduction process by the operation procedure shown in FIG. 13. In FIG. 
14, a polygon model M 0 is an original polygon model. This polygon model M 0 is formed by arraying 28 vertices v v 
v 2 , ...Vgg in 7 rows such that the numbers of the vertices in the respective rows are (3, 4, 5, 4, 5, 4, 3) and connecting 
adjacent vertices. 

[0090] Edges e1 , e2, e3, e4, e5 are selected as such to be successively shrunk by the operation procedure shown 
15 in FIG. 13 to simplify the description since they are not mutually influenced in the case that the respective edges ei (i 
= 1,2, ...5) are shrunk. 

[0091] Accordingly, the data reduction process of FIG. 1 4 shows a case where the edges e1 , e2, ...e5 of the polygon 
model M 0 are successively shrunk by the operation procedure of FIG. 13. Polygon models M.,, M 2 , M 3 , M 4 , M 5 are 
approximated polygon models obtained when the polygon model M 0 is successively approximated by successively 

20 shrinking the edges e1 , e2, ...e5. 

[0092] According to the procedure of the data reduction processing shown in FIG. 1 3, when the polygon model M 0 
is first approximated to the polygon model M 1 by shrinking the edge el by the first edge shrinking processing, edges 
formed by a converged vertex v a and its neighboring vertices v 9 , v 10 , v 1v v 13 , v 1Sl v 18 , v 19 , v 20 (vertices located within 
a dotted-line frame in the polygon model Mj) are influenced. Thus, in the next edge shrinking processing, the estimation 

25 values e q are sorted again for the edges (those located outside the dotted-line frame of the polygon model M 1 ) excluding 
the above edges, and the polygon model M 1 is approximated to the polygon model M 2 by shrinking the edge e2 cor- 
responding to the minimum estimation value Minfe^. 

[0093] Further, since the edges defined by a converged vertex V b and the vertices v 1( v 2 , v 4 , v 6 , v 8 , v 10 , v 13 , v^ 
(vertices located within a dotted-line frame in the polygon model M 2 ) are influenced in the next edge shrinking process- 
30 ing, the estimation values e s are sorted again for the edges (those located outside the dotted-line frame of the polygon 
model M 2 ) excluding the above edges, and the polygon model M 2 is approximated to the polygon model M 3 by shrinking 
the edge e3 corresponding to the minimum estimation value MinleJ. 

[0094] The edge shrinking is applied to the edges e4, e5 by the same method, and the original polygon model M 0 
is approximated to the polygon model M 5 and the number of data is reduced by 5 by the five edge shrinking processings. 

35 [0095] If a process for approximating the original polygon model M 0 to a 6 th polygon model M 6 (not shown) is con- 
sidered, no estimation value is recalculated for the edges having been influenced by the edge shrinking according to 
the operation procedure shown in FIG. 1 3. Thus, in the example of FIG. 1 4, when the estimation value is first calculated 
for 66 edges forming the original polygon model M 0 , the evaluation values are recalculated for 50 edges forming the 
polygon model M 5 without calculating the estimation values again for the approximated polygon models M 1 to M 4 . 

40 Therefore, the estimation values are recalculated 50 times until the & h polygon model M 6 is obtained. 

[0096] On the other hand, since the estimation values are recalculating for the edges influenced by the edge shrinking 
according to the operation procedure shown in FIG. 7, the estimation values are recalculated for 1 6 edges influenced 
by the edge shrinking for each of the approximated polygon models M 1 to M 5 . Thus, according to the operation pro- 
cedure shown in FIG. 7, the estimation values are recalculated 80 (16x5) times until the 6 th polygon model M 6 is 

45 obtained. Therefore, the operation procedure shown in FIG. 1 3 has an advantage of high-speed processing since the 
estimation values are recalculated a reduced number of times. 

[0097] FIG. 15 shows a modification of the flowchart of FIG. 13. In the flowchart shown in FIG. 15, the edges influ- 
enced by the edge shrinking are made not to be shrunk in the next edge shrinking by replacing the estimation values 
thereof by a predetermined estimation value larger than the permissible value. 

so [0098] Specifically, the flowchart of FIG. 1 5 differs from that of FIG. 1 3 in that Step #1 3a is changed to Step #1 3b in 
which the estimation values of the edges influenced by the edge shrinking are set at e^, and the discrimination criteria 
in Step #15 and the objects in Step #17 whose estimation values are recalculated are changed to "edges whose 
estimation values are 6^", and Steps #1 9 and #21 are added after Step #1 7a to shrink the edges whose estimation 
values are recalculated. Accordingly, only part of the flowchart of FIG. 15 different from that of FIG. 13 is described in 

55 detail below. 

[0099] When the edge corresponding to the minimum estimation value Min{ej} is shrunk in the edge/surface shrinking 
section 323 (Step #11), the estimation values of the edges influenced by this edge shrinking processing are set at the 
predetermined estimation value larger than the permissible value (Step #13b) and then Step #7 follows. This 
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operation is performed so as not to shrink these edges in the next shrinking processing, i.e., to apply the edge shrinking 
only to the edges not influenced by the . mlnim um estimation value MhM is 

Se!^ 

mm Subsequently the estimation values of the edges influenced by the edge shrinking processing this time are 
1 t VI So #S and then Step #7 follows again. Thereafter, by repeating the above operations, the next edge 

nrocess at e_ every time the edge shrinking processing is performed (loop of Steps #7, #9, #11 , M). 
EST OMhe^her hand, if the minimum estimation value Mlnfe) becomes smalierthan the permissible vdue . p» 
n Step #9) Zte the edge shrinking processing is repeated, it is discriminated whether there .s any edge having the 
^SSSiSXZ (Step #1 5a" S there are such edges (YES in Step #15a), the estimation values e, (s = 1 , 2, ...) 

[f the minimum estimation value Min{eJ is equal to or below than the perm-ssible value (YES in Step #21), Step 
follows in which the edge corresponding to the minimum estimation value MinfeJ » ^runk 
rnioai on the other hand if there is no more edge having the estimation value of ^ (NO in Step #15a) or the 
Km es*^^^^ the estimation values calculated forthe edges having the estimation value of 

^™eTSger 7m the permissible value (NO in Step #21 ) while the edge shrinking processing is repeated, the 

FIG 15 as in the operation procedure of FIG. 13 by setting the estimation values of such edges at the estimation value 
H 2^^ST?Zh. value, the estimation values are recalculated a reduced number of times and high- 

^oll^^u^Zdl reduction processing shown in F.GS. 13 and 15. the estimation value, are not 
^HJ!^7!^S^ by the edge shrinking processing, i.e.. these edges are not shrunk untol toe edge 
^^DreawiSb prohibited for all edges. However, the prohibition of the edge shrinking processing may be 
leased 9 «7!^cSSSXL such as when the edge shrinking processing is repeated a specified number of 
fn^ofwhe i a difference between the minimum estimation value and the permissible value hes with n j a specified 
ange. and the edge shrinking may be made possible again by recasting the ^^TSSlSSS^ 
raiori Such an arranqement enables the data reduction processing to be applied to the edges for which the edge 

ThorofnrA data can be maximally efficiently reduced within a permissible range. 

^Sm^S^£S!S!Svm data reduction processings shown in FIGS. 1 3 and 15 are designed to increase 
EoceSng speed by mducing the number of times the estimation values are recalculated, the estimation values 

SSr SSil 'values are changed by the edge shrinking or ^^^^ 
to be recalculated after shrinking, for example, those disclosed in the aforementioned paper Surface Simplification 

— 1 - ° btained ? the rsEss 2: 

measuring apparatus 1 is described in the foregoing embodiment, the present invention is not United thereto and is 
widely applicable to three-dimensional data created by computer graphics or the 

roiioi As described above the three-dimensional shape data are reduced by converging two or more vertices to 
ona Ite^t shrink ^the edges or surfaces. Specified estimation values are calculated for all the surfaces defended 
bv shrink th ^Ses or surtais of the polygon model based on the distances between the respective surfaces and 
2 the oSaf vertices LoSn the deformation of these surfaces. The three-dimensional shape data are reduced 
by drinking when the estimation values are equal to or below a 

yalue ASorTgiy the error can be controlled by the estimation standard intuitfceiy easily understandable, i.e.. by 
Snto 2 zig^S ^suTace ir I the polygon model, the data reduction processing can be suitably performed whiiethe errors 

S?iT"~^ 

easilv calculated only by storing an error matrix calculated for the original three-d.mensional shape data. Thus, a high 
^Z Xe^VpL^e and a memory necessary for the calculation of the data reduction processing can be 
efficiently utilized. 
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[0113] Furthermore, the three-dimensional shape data is reduced by repeating the edge or surface based on esti- 
mation values. However, the application of the shrinking to the edges influenced by the previous edge shrinking is 
prohibited. Accordingly, recalculation of estimation values can be executed a reduced number of times, which thus 
enables a higher speed processing. 

[0114] The converged vertex portions, which are influenced by the edge or surface shrinking, are given with a pre- 
determined value to place these portions out of the shrinking, which can thus assure easier shrinking prohibition to a 
particular portion. 

[0115] Moreover, in the case that a specified condition is satisfied, the prohibition is released to allow the shrinking 
processing to the prohibited portion again. Accordingly, data reduction can be accomplished as efficiently as possible 
within a permissible range. 

[01 16] As described above, a method for reducing three-dimensional shape data, comprises the steps of: calculating 
estimation values for portions to be deformed by converging two or more vertices of a polygon model; reducing the 
number of data for the polygon model by converging two or more vertices of the polygon model when the estimation 
values are equal to or below a predetermined permissible value. 

[0117] The portions to be deformed may be surfaces of the polygon model, and estimation values for the surfaces 
are calculated based on distances between the respective surfaces and all the original vertices involved in the surface 
deformation. 

[0118] The estimation value may be an error c defined by Equation (A): 



P f denotes a column vector for an equation of a plane Dj subject to the deformation; 
F{.denotes an inverse matrix of P t ; 

KVy denotes an error matrix for a j-th vertex Vj involved in the deformation of the plane p, when an error matrix for 
data of an original vertex v involved in the deformation of the plane p, is defined by Kv =V- ^(V denotes a column 
vectortfor.f the data); 

m, denotes a sum of error matrix Kvy for a vertex Vj involved in only the deformation of the i-th plane p,; and 
n denotes a sum of planes p, to be deformed. 

[0119] The polygon model may include a number of triangular polygons. 

[0120] The number of data for the polygon may be reduced by converging two or more vertices of one portion of the 
polygon model after another portion, repeatedly. Before each data reduction, the portion having been involved in all 
the previous data reductions may be defined as a reduction prohibition area. A succeeding data reduction may be 
applied to a portion other than the reduction prohibition area. 

[0121] The reduction prohibition area may be released if a predetermined condition is satisfied. 

[0122] The predetermined condition may be a state that there is no portion to be deformed by converging two or 

more vertices. 

[01 23] The predetermined condition may be a state that a predetermined number of data reductions are completed. 
[0124] An apparatus for reducing three-dimensional shape data, comprises: estimation value calculating means for 
calculating estimation values for portions to be deformed by converging two or more vertices of a polygon model; and 
data reducing means for reducing the number of data for the polygon model by converging two or more vertices of the 
polygon model when the estimation values are equal to or below the predetermined permissible value. 
[01 25] The estimation value calculating means may calculate estimation values for surfaces to be deformed by shrink- 
ing edges or surfaces based on distances between the respective surfaces and all the original vertices involved in the 
surface deformation. 

[0126] The estimation value may be an error e defined by Equation (A): 



n 



m, 




wherein: 
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...(A) 



I"* 

wherein: 

P, denotes a column vector for an equation of a plane p, subject to the deformation; 

P.t denotes an inverse matrix of P,\ H«f„,m«ti«n of the Diane d, whan an error matrix for 

KV, denotes an error matrix for a j-th vertex v, involved In the vdenot es a column 

data of an original vertex v involved in the deformation of the plane p, Is defined by Kv- V V ( v 

"Ten^stm of error m atrix Kv y for a vertex Vj involved in on^me delation of the i-th plane 

n denotes a sum of planes p ( to be deformed. 

[0127] Theestimation^^ 

tuting an original polygon model. When two or more verhces v ^^"^XvJSned Equation (A) by using 
estimation value catenating m 

a sum of error matrix Kv p for each converged vertex v p as a erro ^0^ diagonal elements and 

Termr e in accordance with the above-defined Equation (A) by usmg the restored error matnx. 
[01 29] The estimation value calculating means may calculate 

I"" 



(B) 



[0130, Theestimationvaiuecaicu^ 

[ion (A) by using a weighted error matrix w f • Kv,for the vertex V, m place of Kv, 

[0131] The polygon mode, may include a number 0 polygon model by converging two or more 
[0132] The data reducing means may reduce the number of ^ tor tne p yg r 
tertices of one portion of the polygon mode, after an°^r P ort » n ; ?^on of the p^on model is permissible for 
means for iudging based on a <*^"?«^ een .nvoived in all the previous 

p3, aPP The d prohibrtion area defining means may change the estimation value for an involved portion to a predeter- 

s^r^ - - area a ~ ordins to a 

PT^SS ned condition may be a state that there no portion to be deformed by converging two or 
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more vertices. 

[01 36] The predetermined condition may be a state that a predetermined number of data reductions are completed. 
[0137] The judging means may include: minimum estimation value calculating means for calculating a minimum 
estimation value; comparing means for comparing the minimum estimation value with a permissible value; and deter- 
mining means for determining a portion having the minimum estimation value as a portion to be converged when the 
minimum estimation value is equal or below the permissible value. 
[0138] The portion to be converged may be an edge of a triangular polygon. 
[0139] The portion to be converged may be a surface of a triangular polygon. 

Claims 

1 . A method for reducing three-dimensional shape data, comprising the steps of: 

is calculating estimation values for portions to be deformed by converging two or more vertices (v,) of a polygon 

model (Mj); 

reducing the number of data for the polygon model (Mj) by converging two or more vertices (v,) of the polygon 
model (Mj) when the estimation values are equal to or below a predetermined permissible value. 

20 2. A method according to claim 1 , wherein the portions to be deformed are surfaces of the polygon model (Mj), and 
estimation values for the surfaces are calculated based on distances between the respective surfaces and all the 
original vertices (v,) involved in the surface deformation. 



25 



3. A method according to claim 2, wherein the estimation value is an error e defined by Equation (A): 



n ft. 



£= — •..(A) 

30 Jjm, 
wherein: 

35 p f denotes a column vector for an equation of a plane p f subject to the deformation; 

P f T denotes an inverse matrix of Pj> 

Kvj denotes an error matrix for a j-th vertex Vj involved in the deformation of the plane p, when an error matrix 
for data of an original vertex v involved in the deformation of the plane P| is defined by Kv = V • \F ( V denotes 
a column vector for the data); 

40 my denotes a sum of error matrix Kvy for a vertex Vj involved in only the deformation of the i-th plane Pj; and 

n denotes a sum of planes pj to be deformed. 

4. A method according to any one of claim 1 to 3, wherein the polygon model (Mj) includes a number of triangular 
polygons. 

45 

5. A method according to claim 1 , wherein the number of data for the polygon is reduced by converging two or more 
vertices (vi) of one portion of the polygon model (M,) after another portion, repeatedly, and before each data re- 
duction, the portion having been involved in all the previous data reductions is defined as a reduction prohibition 
area, and a succeeding data reduction Is applied to a portion other than the reduction prohibition area. 

50 

6. A method according to claim 5, wherein the reduction prohibition area is released if a predetermined condition is 
satisfied. 

7. A method according to claim 6, where the predetermined condition is a state that there is no portion to be deformed 
55 by converging two or more vertices (V|). 

8. A method according to claim 6, wherein the predetermined condition is a state that a predetermined number of 
data reductions are completed. 
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9. An apparatus for reducing three-dimensional shape data, comprising: 

estimation value calculating means (321) for calculating estimation values for portions to be deformed by 

mte^ 
permissible value. 

surfaces and ail the original vertices (v,) involved in the surface deformation. 
11. An apparatus according to claim 10, wherein the estimation value is an error e defined by Equation (A): 

B ...(A) 



35 



40 



45 



50 



55 



wherein: 

P, denotes a column vector for an equation of a plane p, subject to the deformation; 

fo/data of an original vertex v involved in the deformation of the plane p, .s def .ned by Kv = V ■ \T ( Vdenotes 
l^Z^ZZTlr* Kv, for a vertex v, invoived in on* the deformation of the i-th plane p,; and 
n denotes a sum of planes p| to be deformed. 
,2. An apparMu, according t. any oca of clalma 9 to 11 , whomin o,a aa«ma„on «.,u< 'cowing jmans OKI) c* 
c u ,a,2 an „-* K,,o, -~ -ro~,! Z'SSSS^Z^SXR 
ZSZZ£K£3Z> wZTT^na, Eq u2,»n ,A, » umng a .on, o. a™, ma,* Kv, tor 
each converged vertex v„ as an error matrix Kv q for the vertex v q . 

tion (A) by using the restored error matrix. 
14. An apparatus according to claim 11 , wherein the estimation value calculating means (321) calculates 
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1 5. An apparatus according to claim 1 1 , wherein the estimation value calculating means (321 ) calculates an error e in 
10 accordance with the above-defined Equation (A) by using a weighted error matrix Wj -Kvy for the vertex Vj in place 

of KVy. 

16. An apparatus according to any one of claims 9 to 1 5, wherein the polygon model (M ( ) includes a number of triangular 
polygons. 

15 

17. An apparatus according to claim 9, wherein the data reducing means (323) reduces the number of data for the 
polygon model (Mj) by converging two or more vertices (v { ) of one portion of the polygon model (M,) after another 
portion, repeatedly, further comprising: 

20 judging means forjudging based on a calculated estimated value whether a portion of the polygon model (Mi) 

is permissible for the deformation; and 

prohibition area defining means for defining a portion having been involved in all the previous data reductions 
executed by the data reducing means as a prohibition area to keep a succeeding data reduction from being 
applied. 

25 

18. An apparatus according to claim 17, wherein the prohibition area defining means changes the estimation value 
for an involved portion to a predetermined value to renderthe involved portion impermissible by the judging means. 

< r 19. An apparatus according to claim 1 7 or 1 8, further comprising prohibition releasing means for releasing the prohi- 

30 bitton area according to a predetermined condition. 

20. An apparatus according to claim 1 9, where the predetermined condition is a state that there is no portion to be 
& deformed by converging two or more vertices (v,). 

35 21 . An apparatus according to claim 1 9, wherein the predetermined condition is a state that a predetermined number 
of data reductions are completed. 

22. An apparatus according to any one of claim 1 7 to 21 , wherein the judging means includes: 

40 minimum estimation value calculating means for calculating a minimum estimation value; 

comparing means for comparing the minimum estimation value with a permissible value; and 

determining means for determining a portion having the minimum estimation value as a portion to be converged 

when the minimum estimation value is equal or below the permissible value. 

45 23. An apparatus according to any one of claims 1 7 to 22, wherein the portion to be converged is an edge of a triangular 
polygon. 

24. An apparatus according to any one of claims 17 to 23, wherein the portion to be converged is a surface of a 
triangular polygon. 

50 



55 
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